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Summary 


TABLE I 


An accurate determination has been made of the unper- 
turbed atomic hydrogen ground state hyperfine transition fre- 
quency (F=l,m=0 - F=0 , m=0) in reference to present world 
wide realizations of internationally defined time interval. In 
relation to the international atomic time system as derived by 
the Bureau International De L'Heure, the composite value is 
1,420, 405, 751. 7755±0, 0031 HZ, 


Introduction 


The coherently synthesized integral 1 Mhz output frequency 
of an atomic hydrogen maser (NP-4) developed at NASA's God- 
dard Space Flight Center (GSFC) 1 was compared with the local 
atomic time scale system USNO(MEAN) derived by the United 
States Naval Observatory (USNO) from the mean of selected ce- 
sium standard clocks. 2 The oscillation frequency of the hydro- 
gen maser is related through determinable small correction 
factors to the unperturbed ground state hydrogen atom hyperfine 
transition frequency (F=l , m=0 - F=0 , m=0) which we hereafter 
designate f H . The USNO (MEAN) time scale is compared with 
the local standard frequencies of seven independent national 
time standards laboratories which contribute to international 
atomic time (IAT) as derived by the Bureau International 
De L'Heure (BIH) 3 . The period of the measurement is Febru- 
ary 2, 1972 to April 18, 1972. 


Absolute Hydrogen (F=l,m=0 -F=0,m=0) Frequency with 
Respect to Time Scales of World Time Standard's Laboratories 


F h = 1,420,405,751. + .77XY HZ 


F 

(France) 

USNO 
(Wash. DC) 

IATt 

(BIH) 

NRC* 

(Canada) 

PTB* 

(W. Germany) 

RGO 

(UK) 

ON 

(Swiss) 

NBS 

(Boulder Col.) 

0.7753 

0.7753 

0.7755 

0.7755 

0.7757 

0.7762 

0.7764 

0.7772 


Coordination Precision 0.000,08 Hz 


NP-4 Inaccuracy 0.0024 Hz 


♦Maintain laboratory basic standard - see Appendix IV. 
t Composite system. 


The data and calculations involved in Table I are given in 
the Appendix. These include net absolute corrections to the 
NP-4 oscillation frequency of -0.0879 Hz with the NP-4 associ- 
ated error of ±0.0024 Hz (1 a). The value of f H with respect to 
IAT (BIH), which represents the coordinated international atomic 
time system, is taken as the resultant final value of f H . The 
error assigned to IAT (BIH) is discussed in Appendix IV. The 
final results are as follows : 


f H = 1420,405,751.7755Hz 
a (NP-4) * ±0.0024 Hz 


The opportunity for the present determination arose as the 
result of a cooperative program between the USNO and GSFC to 
evaluate the long term frequency stability of the NP-4 hydrogen 
maser with respect to the USNO time scale system. The USNO- 
(MEAN) scale is based upon the world's largest maintained ce- 
sium clock ensemble. The stability comparison began on Janu- 
ary 28, 1972 when NP-4 was placed at the USNO, and it is 
anticipated that uninterrupted comparisons will continue through 
most of 1972 at the minimum. 

The possibility of making the NP-4 absolute frequency 
measurement with regard to several international standards 
laboratories was realized soon after the stability comparisons 
began when it was observed that the errors associated with 
today's international time coordination links over extended 
measuring periods (two months in the present case) are small 
in relation to inaccuracy factors or reproducibility factors in 
presently used atomic time standards. 

Results 

The principal results are given in Table I, which gives f H 
as related to the NP-4 maser with respect to the standard clock 
systems of the laboratories associated with the BIH. 


a(IAT) = ±0.0020 Hz. 

The net resultant error is taken as a 2 = a 2 (NP-4) + a 2 (IAT) 
whence a =±0.0031 Hz. 

This value for f H is compared with other recent experimen- 
tal determinations in Table n . 


TABLE H 


f H = 1420,405,751. +.7XYZH z 


Laboratory (Date) 

f H - Hz 

Inaccuracy la Hz 

LSRH (69)4 

-0.7782 

0.0036 

NASA (70)5 

-0.7767 

0.0031 
(Appendix V) 

NBS (70)« 

-0.769 

0.002 

NBS (70)6 

-0.767 

0.002 

Morris (70-71) 7 

-0.770 

0.003 

IAT (Present) 
Paper 

-0.7755 

0.0031 



Conclusions 


4.7 Hz RMS. The magnetic field is given by 


The data and results presented point to several significant 
conclusions. 

1. The agreement of NP-4 absolute frequency with results 
of previous determinations is quite good (Table IT). It should be 
noted that no specific experimental effort nor interruption of 
normal real time operation of NP-4 as a clock or frequency 
standard was required to obtain the data. 

2. The present results form the basis for later compari- 
sons of NP-4 frequency and stability with regard to the USNO 
and to world vide standards laboratories, which will be useful 
in evaluating the very long term hydrogen maser performance. 

3. It is possible through present world -wide time coordi- 
nation links (Loran-C, Omega, Satellite transponder, moon 
bounce, VLF, portable clocks, etc.) to make precise time and 
frequency comparisons among videly dispersed international 
standards with a precision which exceeds present atomic stand- 
ard accuracy or reproducibility by well over an order of mag- 
nitude for measuring times as short as two months. Timescale 
coordination on either a local or world wide basis is much more 
easily achievable today than absolute stability or intrinsic ac- 
curacy (Table I). 

Appendix 

I. Oscillation frequency of NP-4 with respect to USNO(MEAN) . 

For the 75 day period Feb. 2, 1972 to April 18, 1972 the 
operational fractional offset of NP-4 (1 MHz) with respect to 
USNO(MEAN) was -10.3 x 10‘ 13 ±0.2xl0' 13 (NP-41ow). The 
error given is the 1 o variance (N = 2, T = 5 Days). There 
was no significant change in relative offset during the measure- 
ment period. The oscillation frequency, f n , of NP-4 is related 
by coherent synthesis to the 1 MHz output f { by the relation 

fn - N,f, 

Where N s = 1,420.405,751,688,8 
thus Af n = N s Afj = N s (-10.3 x 1 0~ 1 3 ) ( 1 0 6 ) 

= -0.001463 Hz 

thus f n = N s [USNO(MEAN) + Afj]. 

Addition gives the oscillation frequency of NP-4 : f n = 1 , 420 , 
405,751.687,34 Hz with respect to USNO (ME AN). 

n. Absolute corrections for NP-4 oscillation frequency. 

The absolute (unperturbed) frequency f H , of the hydrogen 
atom transition frequency is related to the measured value f n 
by the relation 

f n = fji + Afpj + Afp + Af c + Af w . 

Where A f H = Magnetic field correction 

Af D = second order doppler correction 
Af c = Cavity plus spin exchange correction 
Af w = Wall shift correction 

1. Magnetic field correction. 

The Zeeman frequency, f z , of NP-4 was checked weekly. 
The average value was 979 Hz with a weekly deviation of 


H = , A J Z , ^ = (699.3 ± 3.4) x 10' 6 gauss. 

1.40 x 10 6 

The magnetic correction is given by Af H « 2750 H 2 whence 
it is calculated: Af H = (0.001,345 ± 0.000,001,3) Hz. 

2. Second order doppler correction. 

Af D = -3 . 9557 x 10" 4 T where T = absolute temperature of 
the maser bulb (atom temperature). For NP-4, T = 309.1 ± 

0. 2°K, whence it is calculated: Af D = (- 0.060,451 ±0.000, 

039) Hz. 

3. Cavity and spin exchange correction. 

The NP-4 maser is automatic aHy tuned continuously 1 in a 
manner which removes the composite cavity pulling and spin 
exchange correction. The NP-4 cavity correction voltage is re- 
corded continuously. The present correction deviations in addi- 
tion to previous intercomparison data on this design substantiate 
the following estimated error, 

Af c = 0 ± 5 x 10' 14 (Max). 

4. WaU shift correction. 

The wall shift parameters used in calculations of correc- 
tions and errors have been substantiated by several unpublished 
intercomparisons of NASA (GSFC) designs and early Varian 
Associates H-10 designs. All five GSFC H-masers as weU as 
two Varian H-10 masers give consistent frequency agreement 
within 1 x 10* 12 RMS. Small 5-1/2 inch diameter bulbs coated 
at Varian were used originally in the H-lO’s, as well as the 
same masers with bulbs recoated at GSFC. The NASA masers 
as weU as the Varian masers were consistently processed with 
the same materials and techniques as originaUy used in experi- 
ments with wall shift correction at Varian Associates®* The 
wall shift equation is given by 

Af W “ f V 11 - a ( t ‘ 4 °)] 

Where A/V is the surface to volume ratio for the maser storage 
bulb, K = 0.208 Hz-in, a = 5 x 10' 3 and t is the bulb tempera- 
ture in degrees Celsius. For NP-4 the bulb is a 5. 549 liter 
cylinder with hemispherical ends. The inside diameter is 
14. 14cm, and the length of straight section is .25. 8 cm. The 
surface to volume ratio is calculated from these figures to be: 

A/V = 0.8147 in* 1 . 

Using t = 35.9°C and the above figures it is found: 

Af w = (-0.02883 ± 0.00241) Hz. 

The net error given is the root sum square of 1 x 10‘ 12 , 
the statistical deviation among individual masers, plus 1. 2 x 
10" 12 * , the estimated maximum error in the original wall 
shift experiments. This results inl,6xl0" 12 assigned wall 
shift uncertainty. 


*The original published uncertainty was 0.43 x 1 0” 1 2 , but this has been increased 
to the stated value upon the present authors judgment of originally unanticipated 
experimental uncertainties. 



5. Net maser absolute correction and absolute frequency 
f H vs USNO(MEAN). 

Summing the corrections in 1 through 4 gives 
f n = + Af (Total) 

Af (Total = -0.08794 ±0.00241. 

Here the net error is taken as the root sum square of the 
individual contributions, wherein the wall shift correction is 
the preponderant effect. 

Subtracting Af (Total) from f n gives, with respect to USNO- 
(MEAN), f H =1,420,405,751. 775, 28 Hz with ±0.00241 Hz un- 
certainty assigned to NP-^1 alone. The non- significant last 
digit is retained for later computational purposes only. 

HI. USNO(MEAN) and f H in relation to other independent local 
time scales. 

Table HI gives time scale information between laboratories 
and the resultant relations with the USNO and NP-4. Interlab- 
oratory data was taken from BIH circulars D65 (April 5, 1972) 
Section 4, and D66 (May 4, 1972) Section 4. 


The error associated with an individual time coordination 
with IAT is conservatively estimated from the BIH circular data 
as 0. 2 microseconds. Thus the net error becomes \J~2 larger 
with both end periods taken into account, or approximately 0.3 
microseconds. The fractional individual coordination error is 
thus found to be ±5.8 x 10* 14 , or ±0.000, 082 Hz at the hydrogen 
frequency. To obtain f H , the absolute frequency, in relation to 
the individual time scales, the offset between USNO and IAT in 
Table HI was added to Af H in the fifth data column, and these 
numbers were then added algebraically to f H (USNO) ME AN. 


IV. IAT inaccuracy assessment. 

Only two national laboratories are presently using basic 
laboratory long beam cesium standards with published accuracy 
capability in the present range of interest to help generate their 
time scales. These are NRC (Canada), 10 with an error limit 
of 1 1.5 x 10~ 12 , and PTB (W. Germany) 11 with a stated error 
limit of ±4,1 x 10" 13 . The U. S. National Bureau of Standards 
began a major modification program in 1969 to improve the NBS 
3. 66 meter primary frequency standard NBS HI, but the new 
standard, which has been designated NBS V, is not yet opera- 
tional, and the NBS time scale accuracy is thus specified at 
present as 1 x 10* 12 , of which 5 x 10~ 13 is allocated to the last 


TABLE TR 


Laboratory 

IAT 

(Microseconds) 
Feb. 2, ’72 

April 2, '72 

Change 
p SEC 

Fractional 
Parts 
in 10 13 * 

Equivalent t 
Af H -MilliHz 

Af H WRT 
USNO (MEAN) 
-Mill Hz 

f H = 1420,405, 
+0. 77XYZ Hz 

F (1) 
(France) 

- 56.2 

- 57.0 

-0.8 

- 1.54 

-0.219 

-0.027 

0.77525 

USNO (2) 
(Wash. D.C.) 

-34398.5 

-34399.2 

-0.7 

- 1.35 

-0.192 

0 

0.77528 

IAT (3) 
(BIH) 

0 

0 

0 

0 

0 

+0.192 

0.77547 

NRC (4) 
(Canada) 

1.2 

1.1 

+0.1 

+ 0.19 

+0.027 

+0.219 

0.77550 

PTB (5) 

(W. German) 

- 377.4 

- 376.5 

±0.9 

+ 1.74 

+0.247 

40.439 

0.77572 

RGO (6) 
(England) 

1.7 

1.1 

+2.8 

+ 5.40 

+0.767 

40.959 

1 

0.77624 

ON (7) 
(Swiss.) 

+ 16.5 

+ 19.8 

+3.3 

+ 6.37 

+0.905 

+1.097 

0.77638 

NBS (8) 
(Boulder Col.) 

-45200.9 

-45194.3 

+6.6 

+12.73 

+1.807 

+1.910 

0.77719 


(1) Observatory of Paris, cs clock ensemble, with several cesium clocks in other locations in France, all commercial H. P. , or 
with H. P. tubes. 

(2) USNO ensemble, 14 cs clocks this measuring period, all commercial H.P. 

(3) Derived Time Scale: weights USNO 25%, NBS 12.5%, NRC 12.5%, F 12.5%, PTB 12.5%, RGO 12.5%, ON 12.5%. 

(4) NRC Ottawa, cs clock ensemble, all commercial H. P. , & 1 basic laboratory standard. 

(5) PTB Braunschweig, cs clock ensemble, all commercial H. P. , & 1 basic laboratory standard. 

(6) RGO Herstmonceux, cs clock ensemble, all commercial H.P. 

(7) Observatory Neuchatel, cs clock ensemble, all commercial H. P. 

(8) NBS Boulder, cs clock ensemble, with several cs clocks in other locations, all commercial H. P. 

(*) The net elapsed time, Feb. 2 to April 2 is 60 days, or 5.184 x 10 6 seconds. 

(f) Fractional change multiplied by f H = 1.420 x 10 9 . 



(1969) inaccuracy assessment of NBS m 13 , and the remainder 
to the long term instability of the NBS commercial cesium clock 
ensemble over the 3 year period since NBS m was operational 13 . 

The other time scale systems relate inaccuracy primarily 
to the reproducibility specifications of the commercial cesium 
standard clocks used, which in turn have been related through 
manufacturers inter comparisons with house standards, or with 
occasional further inter comparisons with laboratory long beam 
cesium standards maintained by national laboratories. 

In view of this complex system, the inaccuracy associated 
with the IAT System is arrived at as follows. First, we note 
that the difference between USNO and NBS is 0. 0019Hz (13.4 x 
10‘ 13 fractionally). 

The smallest inaccuracy factor for laboratory cesium beam 
standards officially reported to date is that of PTB (W. Ger- 
many), namely ±4.1 x 10' 13 . 

Thus the net IAT inaccuracy factor is taken as the root sum 
square of 13.4 x 10 3 and 4.1 x 10 " l 3 which gives 

a (IAT) = 14.0 x 10' 13 factionally or 

o (IAT) = 0.0020Hz at the hydrogen frequency. 

The above procedure, while somewhat arbitrary, does not 
seem unreasonable, and is judged to give a realistic inaccuracy 
factor for the overall IAT System. Combining the error above 
with the NP-4 associated error as the root sum square gives 
the net error of a (Net) = ±0.0031 Hz or ± 2. 2 x 10 _1 2 
fractionally. 


V. Correction of NASA (70) Data. 

Reference (5) gives the net error in the f H determination 
using the NASA (NX) hydrogen maser and five HP cesium clocks 
as a = 0.0024Hz. This is judged to be optimistic. The un- 
certainty for wall shift given in Reference (6) for the NASA (NX) 
maser did not include all the uncertainty assumed herein for 
NP-4 wall shift correction. To put the (NX) number on the 
same conservative basis as the present data, the (NX) data 
should be corrected to give for (NX) alone o (NX) = ±0.0024 Hz. 
This may be combined with the net Cs inaccuracy, assumed in 
Reference (5) to be ±1. 33 x 10" 1 2 or a cs = ±0. 0019Hz. Taking 
the root sum square of these gives o (Net) = ±0.0031 Hz. This 
is the value the present authors assume in Table II . 
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